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We study the dynamical structure factors of the spin-1/2 XXX chain at finite magnetic field h, 
focusing in particular on the singularities at excitation energies in the vicinity of the lower thresholds 
of the leading-order dominant excitations. We derive the exact momentum and spin-density 
dependences of the exponents C^{k) controlling the singularities for both the longitudinal (r = 1) 
and transversal (r = t) dynamical structure factors for the whole momentum range k £ [0, tt] , in the 
thermodynamic limit. In that limit we argue that the higher-order excitations change neither (4^{k) 
nor the sharpness of the singularities. We discuss the relation to neutron scattering and suggest 
new experiments on spin-chain compounds using a carefully oriented crystal. 

PACS numbers: 75.10.Pq, 71.10.Pm 


The anisotropic spin-1/2 Heisenberg XXZ chain in a 
magnetic field h with anisotropy parameter A > 0 and 
exchange integral J is a paradigmatic example of an iii- 
tegrable strongly correlated system [^. Its Hamiltonian 
is, 

^ ( E ^) +9FBhj2Sh 

j=l \T=x,y / j=l 

( 1 ) 

where for simplicity we take g = 2, gs is the Bohr mag¬ 
neton, and are components of the spin-1/2 oper¬ 
ators at site j = We consider states with spin 

S s]0, L/2] and spin density m = 2S^jL S]0,1] where 
= ^{L — 2M) is the spin projection and M the num¬ 
ber of down spins. 

At the isotropic ^int, A = 1, the XXZ chain becomes 
the XXX chain |2|,[3|| and contains antiferromagnetic cor¬ 
relations that have been observed in dynamical quantities 
measured in experiments on spin-chain compounds 
For instance, its dynamical structure factors. 




dte-“‘(0|A“(t)5“(0)|0) 


( 2 ) 

where a = x,y, z are important examples of relevant 
quantities from the experimental viewpoint [^. 

The isotropic A = 1 model also poses the most chal¬ 
lenging technical problems for theory. This holds true 
also for h = 0, which is the case that has been most 
heavily investigated im. Indeed, as shown in Ref. 
[ll| . for A G [0, 0.8] and h = 0 nearly all the S^^{k,oj) = 
S^^{k, uj) integrated intensity spectral weight stems from 
two-hole Bethe ansatz (BA) excitations whereas as A in¬ 
creases from 0.8 to 1.0, the four-hole excitations increas¬ 
ingly contribute as A —^ 1. Specifically, in that limit the 
contributions to S^^(k,w) from the S' = 1 and S^ = 0 


two-hole excitations correspond to a relative integrated 
intensity of k, 0.75 Q. This increases to k, 0.99 if one 
considers in addition S = 1 and S^ = 0 four-hole ex¬ 
citations. Importantly, however, all the S^^{k,uj) sin¬ 
gularities are determined by contributions from two-hole 
excitations and the exponent that controls them can be 
shown from purely phenomenological considerations to 
be fixed at h = 0 by the SC/(2) invariance alone 0. 


At h > 0 one has that S^^{k,w) ^ S^^{k,w) = 
Syy{k,oj). Previous studies on these dynamical struc¬ 
ture factors in the XXX model refer to finite systems 
and rely on numerical diagonalizations |l2|, e valuation 
of matrix elements between BA states [Sl, [iSl, [14 j, and 
the form-factor method [l^, which applies to the XX Z 
chain 00. The studies of Refs. [ 13 , [0 refer mostly 
to A < 1. The line shape in the vicinity of the dynamical 
structure factor thresholds is predicted to be controlled 
by momentum dependent exponents (l8| . but that de¬ 
pendence remains an important open problem. 


On the other hand, our present study and results refer 
to the thermodynamic limit (TL). Specifically, we clar¬ 
ify one of the unsettled questions concerning the physics 
of the spin-1/2 XXX Heisenberg model at finite mag¬ 
netic field h by deriving the exact momentum and spin- 
density dependences of the exponents that control the 
lower thresholds singularities in S^^{k,uj) and S^^{k,uj). 
We discuss the relation of our theoretical results to the 
w dependence of the magnetic scattering intensity ob¬ 
served in inelastic neutron scattering and suggest new 
experiments using a carefully oriented crystal. 

Our calculation approach uses the pseudofermion dy¬ 
namical theory (PDT) developed in Refs. 0 0 for 
the one-dimensional (ID) Hubbard model in the TL [0, 
which can also be used for the XXX chain. When 
the lower thresholds of the spectral and dynamical func¬ 
tions refer to the branch lines as defined within that the- 
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ory, the exact line shape in the vicinity of such thresh¬ 
olds is indeed controlled by momentum-, density-, and 
interaction-dependent exponents. These are expressed in 
terms of dressed phase shifts that are known from the 
BA solution. 

After the PDT was introduced, a set of novel methods 
were developed to tackle also the finite-energy physics 
of integrable and non-integrable ID correlated quantum 
problems, beyond the low-energy limit 22-23|- For in¬ 
stance, the same exact momentum dependence derived 
by the PDT for the exponents that control threshold sin¬ 
gularities of the ID Hubbard model one-electron spectral 
functions [l^ has recently been obtained in the frame¬ 
work of a mobile impurity approach [22 1. 

For h > 0 the singularities that dominate the line 
shape for small excitation energy (w — uj'^{k)) near the 
lower thresholds w£(fc) of the longitudinal (r = 1) and 
transversal (t = t) dynamical structure factors S^^{k,uj) 
and S^^{k,uj), respectively, Eq. ([2]), are within the PDT 
determined by class (ii) excitations generated by specific 
leading-order processes in the BA distributions. (Class 
(ii) excitations are |S'^| = S excited states Q.) At the 
present isotropic A = 1 point, higher-order class (ii) ex¬ 
citations generated by additional higher-order processes 
in the BA distributions also contribute to the dynamical 
structure factors but in the TL they change neither the 
momentum dependent exponents nor the corresponding 
singularities sharpness. Moreover, class (ii) excitations 
described by complex BA rapidities are gapped for h > 0 
and except for very small h have nearly vanishing spec¬ 
tral weight. For instance, at m = 0.5 their contributions 
correspond to a relative intensity of about 3 x 10~^ for 
anisotropy A = 0.3 and 4 x 10“^ for A = 0.7 [l3|. Even 
for A = 1, the corresponding estimated relative intensity 
is not larger than 10“®. Our study focuses mainly on the 
m > 0.15 range for which their contribution is negligible. 

The lower threshold singularities of the dynamical 
structure factors, Eq. are for h > 0 determined by 
excitations whose BA rapidities Xj are real. There is a 
one-to-one correspondence, Xj = X{qj), to the BA quan¬ 
tum numbers qj = ^Ij where j = 1, • ‘ ‘ ,L — M and Ij 
are successive half-odd integers for L — M even and inte¬ 
gers for L — M odd. Out of the set of j = 1, • ‘ ‘ ,L — M 
quantum numbers {qj}-, a subset of a = 1, • • • , M num¬ 
bers {qq,} are occupied. These refer to the "BA band 
particles" whereas the remaining L — 2M = Lm values 
qj correspond to the "BA band holes", which here are 
simply called particles and holes, respectively. The mo¬ 
mentum of these states is k = M tt + 
numbers {qj} indeed play the role of BA band momen¬ 
tum values. 

We consider the TL, L —^ oo, within which the set 
oi j = 1, • • • ,L — M momentum values {qj} (such that 


qj+i — qj = ^) may be replaced by a continuum mo¬ 
mentum variable, q € [—^Ft]? ^tnd the rapidities 
Xj = X{qj) by a rapidity function, A = X{q) € [— 00 , 00 ], 



with A(±A:Ft) = ± 00 . Except for 0(1/L) corrections, 
one has that = §(1 -I- m) and kp^ = ^(1 — m) with 
q G [—kpi, kpi] corresponding to the occupied momen¬ 
tum values of TO S [0,1] ground states. Hence for to > 0 
such states are populated by holes for jgj G [kpi,kpj]. 
The BA band is full for the m = 0 ground state. (Its 
excited states holes are usually identified with spin- 1/2 
whereas those of to > 0 states are 
.) We denote by Xo(q) the ground 


spmons 

called here holes [2 
state rapidity function such that Ao(g') = —Ao(— 9 ), 
Xo(kpi) = B, and Xo(kpf) = 00 where H = 00 for to —>■ 0 
and H = 0 for TO —>■ 1. It can be defined in terms of its 
inverse function, q = dX2Tra{X), where the usual 

BA distribution 27rcr(A) is the solution of the integral 
equation 27ra(A) = dX' 'V ' 

The spectra, uj'^{k) = w'"(—fc), of the excited states 
that for h > 0 and spin densities to g]0 , 1 ] control the 
leading order contributions to the dynamical structure 
factors involve the BA band dispersion e{q) = e^{q) — 
£°{kpi) where q G [-kpj,kpj], e°(g) = e°(Ao(( 7 )), and. 


e-°(A) 

$(A,A') 


9 T y 

1 + (2A)2 l + {2X'r 


arctan(A — A') 


'-B 

rB 


dX" 


<J>(A',A) 
4>(A",A') 


TT l-t(A-A")^ 


.(3) 


The latter integral equation defines the rapidity- 
dependent dressed phase shift 27r$(A, A') in units of 2p. 
The group velocity reads v{q) = de{q)/dt = de°{q)/dt, 
with the Fermi velocity v = v{kpi) playing an important 
role in the low-energy physics. The dispersion e^{q) con¬ 
trols the spin density curve, him) = —£^{kp\)/{2^B) 
with TO = 1 — 2kpi. The range to G]0, 1] refers to 
h g] 0 , he] where he = is the critical field at which 

fully polarized ferromagnetism is achieved. 

The processes that generate the excited states that 
control the lower threshold singularities correspond to 
specific values of the right (t = I) and left (i = —I) 
Fermi points particle number deviations SNp^. It is more 
convenient though to use the corresponding number devi¬ 
ations SNp = ^i^Ft and SJp = A 

Within the PDT, the exponents that control the line 
shape in the vicinity of the theory branch lines that coin¬ 
cide with T = I, t lower thresholds involve the functional, 

2A:(g)= (^.^+eo<5J^ + C$(.A:Fi,<z)y • (4) 

Here l = ±1, the momentum q is that of a parti¬ 
cle created {C = I) or annihilated {C = —1) un¬ 
der the transitions to the excited states, 4>(g,q') = 
$(Ao((;), Ao(g')) is a momentum dependent dressed phase 
shift, $(A, A') obeys Eq. ([3|), and the related parame¬ 
ter, ^0 = 1 + ®(fcF 4 ,, kp^) - ^{kpi,-kpj,), increases from 
^0 = l/v^ for TO —>■ 0 to ^0 = 1 as TO —>■ I. Indeed, 
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for m —?> 0 the dressed phase shift $(q, g') has the lim¬ 
iting value ^{i,kp^,q) = (,/( 2 -\/ 2 ) for q ^ tkpi whereas 
Hikp^^ikpi,) = Li3-2V2)/{2V2). 

The longitudinal spectrum, uj\k) = uj\—k), refers to 
excited states that are generated by one particle-hole pro¬ 
cesses and conserve M, 

u:\k) =-£{qi) + £{q 2 )-, k = q 2 - qi . (5) 

On the other hand, the spectrum a;‘(fc) of S^^{k,uj) = 
{k,uj) + S ^{k,uj)] is here expressed as the su¬ 
perposition of the spectra (i) (fc) and (ii) uj ^{k). 
The spectra refer to states that are generated by 

a zero-energy 5M = SNpi = ±1 process at a t Fermi 
point along with an overall BA band shift Sqj = 
which renders it a net zero-momentum process, plus one 
(i) particle-hole and (ii) two-hole process, 

a;+“(fc) = £{q2) - £iqi) ; k = n + q2 - qi , 
uj-+{k) = -e((7i) - e(g2); k = tt - qi - q2 . (6) 

All above spectra refer to spin densities m S]0,1] and mo¬ 
mentum k € [ 0 , 7 r]. In Eqs. ([5]) and (jb]), qi € {—kF],,kF\\ 
and q2 G [kFi,kF^] for the longitudinal spectrum, q2 G 
[—fej’-l-, — fc_F 4 ,] for the H— spectrum, and (72 G [— kp^] 
for the —h spectrum. 

In the m —?► 0 limit, the spectra uj’‘{k), Eq. ([5|), and 
0 ;“' (fc), Eq. ([6|), reduce to their lower thresholds as 

h —>■ 0, which in that limit becomes the lower threshold of 
both UJ ^(fc), Eq. ([5]), and the class (ii) two-hole excita¬ 
tions described by complex rapidities whose gap vanishes 
as h —^ 0. At /i = 0 the spectrum w ''(fc) is also that 
of the = 0 and S' = 1 two-hole excitations of class 
(i), which due to a selection rule Q do not contribute to 
the dynamical structure factors at h > 0. Hence upon 
smoothly turning off h there is a large weight transfer 
from |S^| = S class (ii) excitations for h —?> 0 to degener¬ 
ate = 0 and S = 1 class (i) two-hole states at h = 0. 

A particle (and hole) branch line is for the PDT gen¬ 
erated by excitations where one particle is created (and 
annihilated) outside the q = ikpi Fermi points and all 
remaining processes occur at such points. For both spin 
densities m —>■ 0 and m > m* ~ 0.15, the lower threshold 
of w*(fc) (and uj*'{k)) coincides with a hole branch line for 
k G [0, 2fci7’4^] (and fc G [tt — 2kF]^,'K\) and with a particle 
branch line for k G [2kFi, tt] (and k G [0, tt — 2kFi])- For 
0 < m < m* ~ 0.15, the branch-line spectrum uj^{k) (and 
uj^^k)) does not coincide with the hole branch line for a 
small momentum width near fc = 0 (and k = tt.) For 
simplicity, we consider mostly spin densities m —>■ 0 and 
m > m* « 0.15 for which w£(fc) coincides with branch 
lines and the PDT gives the exact momentum and spin 
density dependence of the exponents that control the line 
shape in its vicinity. (Even for 0 < m < m* Ri 0.15 they 
are a good approximation.) 

Interestingly, the use of the PDT reveals that the lower 
threshold singularities of ik,ijj) are those of {k,uj) 


near the particle branch line and of S ^{k,uj) near the 
hole branch line. Accounting for £{±kFi) = 0, the longi¬ 
tudinal S^^{k,uj) and transversal S^^{k,uj) hole branch 
lines spectra are given by 

^hik) = -£{q), T = l,t, 

k = kpi - g G [0, 2kFi], t = 1, 
k = Tr-kFi-qG[Tr-2kFi,TT], T = t, (7) 

where q G [—kFi-,kFi\- The corresponding particle 
branch lines spectra are given by, 

ujp{k) = £{q) , T = l,t, 

k = kpi + qG [2kFi, tt] , t = 1, 

k = TT-kFi + qG[0,TT-2kFi], T = t, (8) 

with q G [kpi, kp^] and q G [—kp^, —kpi] for the I and t 
particle branch lines, respectively. 

For the (7 = 1 particle and (7 = — 1 hole I branch lines 
the Fermi points number deviations to be used in Eq. ® 
are bNp = —C and SJp = For the t branch lines such 
deviations are given by SNp = 0 and SJp = One then 
finds, 

2A((g) = + , 

2A‘(g) = + , (9) 

where q G [—kpi, kpi] for (7 = —1 and r = l,t, q G 
[kpi, for (7 = 1 and t = I, and q G [-kp-f, —kpi] for 
(7=1 and t = t. 

For small positive values (oj — oj'[(k)) in the vicinity of 
ujp{k) > 0, the dynamical structure factors are according 
to the PDT of the form, 

5““(fc,w) = CI(a;-a;E(A:))«^«, G [0, ^], 
e{k) = -l+Y.2A:{q). (10) 

t —±1 

Here a = z ioi t = I, a = x for t = t, (7^ is a k and ui 
independent constant, and the k values that correspond 
to q in 2A[ (g) are those of Eqs. (O and ([5]). 

The exponent ^'^{k) given in Eq. (flU)) does not ap¬ 
ply near the w = 0 lower threshold soft modes such as 
(fcj, 0) where /cg = 2kFi kg = tt — 2kFi (k^ ^)" 

plane. In this case the PDT reaches the same results 
as conformal field theory (CFT). Indeed, near them the 
functionals, Eq. become the conformal dimensions 
of the t = ±1 fields 2A[ = (t(5Ai./(2^o) + 

Here SNp = 0, SJp = 1 for r = I and SNp = ±1, 
SJp = —1 for T = t. While in this Letter we are mostly 
interested in the finite-energy behavior of S°“^{k,uj), for 
completeness we provide its general form near (fcj,0), 
which can be obtained from CFT. Specifically, one finds 
5““(fc,a;) = Cl, {uj)<i 1 near the branch lines for small 
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Figure 1: Two upper panels of each sub-figure (a) and (b): 
The spectra uj\k) and u!*{k) for (a) m = 0.16 and (b) 
m = 0.75. Lower panels: The exponents ^\k) and 
Eq. (IIOII . that control the singularities in the vicinity of the 
lower thresholds of the spectra plotted here as a function of 
k e [0, tt]. 


CO ~ ±v{k — fcj) and = CJ (uj)‘’0 away from 

such lines for small oj ^ ±u (fc — fcj). Here and Cq 
are k and lj independent constants and the exponents 
read C±i = -1 + 2A^ and CJ = -2 -f X;.=±i 2A[. 

In the present L ^ oo limit, the h > 0 higher-order 
particle-hole excitations do not change the exponent and 
corresponding sharpness of the singularities, Eq. 

They may however change slightly the values of the k and 
w independent constants CJ. For finite systems they give 
rise to a small tail for small negative values (w — ujJ (k)). 
Such finite-L effects are seen in Fig. 1 of Ref. [l7| in 
S^^{k,uj) for anisotropy A = 0.7 and momentum values 
k = 37r/4 and k = 7r/2. Due to the finiteness of \ogL/L, 
the peaks appear sharper if one does not account for the 
contributions from the two particle-hole excitations. Our 


results refer to L —>■ oo so that the latter contributions 
do not affect the sharpness of the peaks, Eq. , since 
\ogL/L —^ 0 as L —>• oo. 

The spectrum ijj\k), Eq. ([S]), and the spectrum 
that results from combination of the spectra (fc) and 
w ^{k), Eq. ([6|), along with the corresponding exponents 
and i*'(k) given in Eq. (jlOp are plotted in Fig. 
[T] for spin densities (a) m = 0.16 and (b) m = 0.75. 
The exponent is negative for fc > 0 at any m value 
whereas is negative for a m-dependent range k G 

\kt , tt] . Here kt increases from fct = 0 for m —?> 0 to fct = 
—2 arctan tan(7r/-\/2)) ~ 0.37 tt for m —^ 1. For the 
k ranges for which ^'^{k) < 0, there are lower threshold 
singularity cusps in 5'““(fc,w), Eq. (fTUl) . 

In the m —>■ 0 limit both the t = I, t lower thresholds 
oj'^{k) coincide with the hole branch line for all k values. 
In that limit the corresponding exponents are given 
by = —1/2 for all k values and the lower thresholds 
oj'^ik) coincide with that of the m = 0 two-hole spectrum. 
Consistently, ^'^{k) = —1/2 is also the value of the known 
exponent that controls the line shape in the vicinity of 
the lower threshold of the latter spectrum [s, EEl. 

In the opposite limit, m —>■ 1, the lower thresh¬ 
olds uj'^{k) coincide with the particle branch line for 
all k values and the t = I, t exponents read ^’’{k) = 
2$(0,fc)[l -b $(0,/c)] and (^‘(fc) = -1/2 -b 2$(0, fc - 
7r)[l -b 4>(0, fc — tt)] where the phase shift reads $(0, q) = 
— i arctan tan (I)). In this limit, = 0 (and 

= 1) for fc —>■ 0 decreases to ^^{k) = —1/2 (and 
^*(fc) = —1/2) for /c — >■ TT. The corresponding m —)> 1 be¬ 
haviors refer to a small but finite M/L value. For both 
zero and finite M values reached as h ^ he, the longitu¬ 
dinal spin dynamical structure factor vanishes in the TL 
Q, the superimposed dynamical structure factor being 
for h ^ he dominated by S^^{k,u)). At h = he Eq. (fTOl) 
is not valid, being replaced by the d-function like distri¬ 
bution, S^^{k,u}) = (w — J(1 -b cosfc)) for k G [0,7r]. 

The structure form factors S^^{k,uj) and S^^{k,uj) 
may be investigated separately in > 0 experiments on 
spin-chain compounds by using a carefully oriented crys¬ 
tal. If the crystal is misoriented, or if a micro crystalline 
sample is used, the S^^{k,uj) and S^^{k,uj) spectral fea¬ 
tures should appear superimposed. Such superimposition 
changes the excitations lower thresholds and leads to the 
broadening of the singularities, Eq. (Unj. However, this 
does not occur at = 0, since S^^{k,uj) = S^^{k,uj). 

These two different situations are clearly seen in the 
magnetic scattering intensity measured at zero- and 
finite-field inelastic neutron scattering experiments of 
Ref. respectiveR, on Cu(C 4 H 4 N 2 )(N 03 ) 2 . In Figs. 
2 (a)-(c) of Ref. [3 the theoretically predicted sharp 
cusps at zero-field, S^^{k,uj) = S^^{k,oj) = — 

are clearly seen at different k values. On the 
other hand, the S^^{k,oj) ^ S^^{k,uj) spectral features 
appear superimposed in the finite-field Figs. 2 (d)-(f) of 
that reference, so that only at fc ~ tt is the theoretically 
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predicted sharp cusp clearly visible. 

In summary, we have obtained the exact momentum 
dependence of the exponents that in the TL control the 
line shape singularities in the vicinity of the lower thresh¬ 
olds of the longitudinal and transverse dynamical spin 
structure factors using the PDT [iSliS , in the funda¬ 
mental case of the spin-1/2 XXX Heisenberg chain in a 
field h > 0. We suggest that more demanding h > 0 ex¬ 
periments with a carefully oriented crystal are carried out 
on Cu(C 4 H 4 N 2 )(N 03)2 and other spin-chain compounds, 
thus yielding separately S^^{k,uj) and S^^{k,uj) whose 
magnetic scattering intensities are expected to display 
the singularity cusps theoretically studied in this Letter. 
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the support by the Beijing CSRC and the FEDER 
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respectively. 


[1] M. Takahashi and M. Suzuki, Prog. Theor. Phys. 48, 
2187 (1972). 

[2] M. Takahashi, Prog. Theor. Phys. 46, 401 (1971). 

[3] G. Muller, H. Thomas, H. Beck, and J. C. Bonner, Phys. 
Rev. B 24, 1429 (1981). 

[4] J. Sirker, R. G. Pereira, and I. Affleck, Phys. Rev. B 83, 
035115 (2011). 

[5] M. B. Stone, D. H. Reich, C. Broholm, K. Lefmann, C. 
Rischel, G. P. Landee, and M. M. Turnbull, Phys. Rev. 
Lett. 91, 037205 (2013). 

[6] K. R. Thurber, A. W. Hunt,T. Imai, and F. C. Chou, 
Phys. Rev. Lett. 87, 247202 (2001). 

[7] N. Motoyama, H. Eisaki, and S. Uchida, Phys. Rev. Lett. 
76, 3212 (1996). 

[8] M. Karbach, G. Muller, A. H. Bougourzi, A. Fledderjo- 
hann, and K. H. Mutter, Phys. Rev. B 55, 12510 (1997); 
A. Abada, A. H. Bougourzi, and B. Si-Lakhal, Nucl. 
Phys. B 497, 733 (1997); A. H. Bougourzi, M. Couture, 
and M. Kacir, Phys. Rev. B 54, R12669 (1996). 


[9] J. - S. Caux and R. Hagemans, J. Stat. Mech. P12013 
(2006). 

[10] A. Imambekov and L. I. Glazman, Phys. Rev. Lett. 102, 
126405 (2009). 

[11] J. - S. Caux, H. Konno, M. Sorrell, and R. Weston, Phys. 
Rev. Lett. 106, 217203 (2011). 

[12] K. Lefmann and C. Rischel, Phys. Rev. B 54, 6340 
(1996). 

[13] M. Karbach and G. Muller, Phys. Rev. B 62, 14871 

( 2000 ). 

[14] M. Karbach, D. Biegel, and G. Muller, Phys. Rev. B 66, 
054405 (2002). 

[15] D. Biegel, M. Karbach, and G. Muller, Europhys. Lett. 
59, 882 (2002). 

[16] N. Kitanine, J. M. Maillet, and V. Tetras, Nucl. Phys. B 
554, 647 (1999). 

[17] J.-S. Caux and J. M. Maillet, Phys. Rev. Lett. 95, 077201 
(2005). 

[18] R. G. Pereira, S. R. White, and 1. Affleck, Phys. Rev. 
Lett. 100, 027206 (2008). 

[19] J. M. P. Carmelo, K. Penc, D. Bozi, Nucl. Phys. B 725, 
421 (2005) ; 737, 351, Erratum (2006); J. M. P. Carmelo, 
K. Penc, P. D. Sacramento, M. Sing, and R. Claessen, J. 
Phys.: Cond. Mat. 18, 5191 (2006); J. M. P. Carmelo, 
D. Bozi, and K. Penc, J. Phys.: Cond. Mat. 20, 415103 
(2008). 

[20] J. M. P. Carmelo, L. M. Martelo, K. Penc, Nucl. Phys. 
B 737, 237 (2006). 

[21] F. H. L. Essler, H. Frahm, F. Gohmann, A. Kliimper, 
V. E. Korepin, in The one-dimensional Hubbard model 
(Cambridge University Press, Cambridge, UK, 2005); J. 
M. P. Carmelo, A. H. Castro Neto, and D. K. Camplrell, 
Phys. Rev. B 50, 3667 (1994); J. M. P. Carmelo, A. H. 
Castro Neto, and D. K. Campbell, Phys. Rev. B 50, 3683 
(1994). 

[22] L. Seabra, F. H. L. Essler, F. Pollmann, 1. Schneider, and 
T. Veness, Phys. Rev. B 90, 245127 (2014). 

[23] A. Imambekov, T. L. Schmidt, and L. 1. Glazman, Rev. 
Mod. Phys. 84, 1253 (2012). 

[24] R. G. Pereira, K. Penc, S. R. White, P. D. Sacramento 
J. M. P. Carmelo, Phys. Rev. B 85, 165132 (2012). 

[25] Elsewhere it will be shown that the holes of m > 0 ground 
states and their excited states have scattering properties 
different from those of spin-1/2 objects, such as the m = 
0 spinons. 



